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Abstract. We propose a novel capacity model for complex networks against cascading failure. In this 
model, vertices with both higher loads and larger degrees should be paid more extra capacities, i.e. the 
allocation of extra capacity on vertex i will be proportional to kj , where ki is the degree of vertex i and 
7 > is a free parameter. We have applied this model on Barabasi-Albert network as well as two real 
transportation networks, and found that under the same amount of available resource, this model can 
achieve better network robustness than previous models. 
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1 Introduction 

Many systems in nature and society can be described by 
networks, including biological and social systems, the In- 
ternet, the world-wide web, friendship networks, computer 
networks, metabolic networks, power grids, scientific ci- 
tations, neural networks, and so on. Therefore, complex 
networks have recently attracted considerable attention 
in physics and other fields. Interestingly, many real- world 
networks share a certain number of common topological 
properties, such as small-world and scale-free properties 
1 1 2, 310]. The network robustness is one of the central 
topics in studying of complex networks. Robustness refers 
to the malfunction avoiding ability of a network when a 
fraction of its constituents are damaged. Previous works 
demonstrated that the heterogeneity of a network induces 
a high resilience to random failure while a high sensitivity 
to intentional attacks |5,(3E]- A common failure of many 
networks is cascading failure triggered by the removal of 
vertices or overload breakdown of vertices. The robust- 
ness of complex networks in response to cascading failure 
of intentional attacks has become a topic of recent interest. 
Prior studies have shown that the fault-tolerant capabil- 
ity of network have a great impact on their robustness and 
function [S]lrJ]l7Tl8T^^ 

l22 l lM p 4 p 5 l l26 W . 28, 29, 30, 31 , 32,33,34, 35, 36] . 

Any failures of vertices in general will change the dis- 
tribution of loads. Here the load of a vertex (also called 



betweenness centrality [37 | l38 j ) is defined as the total num- 
ber of shortest paths passing through this vertex. If the 
load at a particular vertex increases and becomes larger 
than its capacity, the corresponding vertex fails. The over- 
load failure leads to redistribution about the loads and, as 
a result, subsequently cascading failures may occur. Be- 
cause of the global redistribution of load, new overload 
failures may be driven by events happening far away. This 
cascading process may stop after a few steps, while it could 
also spread over the entire network. 

The damage caused by cascading failures can be quan- 
tified by the relative size of the largest connected compo- 
nent G, defined as following 



G = N'/N, 



(1) 



where N and AT' are the numbers of vertices in the largest 
component before and after the cascade, respectively. The 
integrity of a network is maintained if G ~ 1, while break- 
down occurs if G « 0. The relative size G also represents 
the robustness of complex networks against cascading fail- 
ures. The cascading failure can be prevented by assigning 
extra capacities to vertices. Since the extending of network 
capacity will bring economic and technique pressure, it is 
important to explore how to rationally allocate the limited 
capacity onto vertices, and efficiently improve the robust- 
ness of network. The capacity of a vertex is the maximum 
load that the vertex can handle. Assuming the capacity 
Gj of vertex i be proportional to its initial load L^. 
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Gj — A, • Li 



(2) 
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Fig. 2. Go and Gmax vs. a. 



Fig. 1. The robustness G vs. 7 under attacking of the highest- 
load vertex on BA network of size N = 5000 and average degree 
(k) = 4 for different a. 



where A^ > 1 is the tolerance parameter. Excess term (A* — 
1) is the extra capacity of vertex i, which simultaneously 
reflects i's ability of tolerating the additional burden and 
the extra cost to protect i. 



2 Model 

Among the previous works, two models should be paid 
special attention: Motter-Lai (ML) model [10] and Wang- 
Kim (WK) model [33] . ML model assumes the capacity 
Ci of vertex i be proportional to the initial load Li as 



d = (l + a) -Li, i= 1,2, 



(3) 



where a > is the control parameter representing the 
extra capacity. In WK model, the capacity Ci is 

d = (1 + a ■ 0(-fi- - /3)) ■ L h i = 1, 2, . . . , N, (4) 

where &{x) = 0(1) for x < 0(> 0) is a two- valued function, 
namely the Heaviside step function, Lmax = niax,(L,-), 
a G [0, 00) and f3 £ [0, 1] are two control parameters. Here, 
each vertex can be in one of two states: assigned some ex- 
tra capacity, or not. When (3 — 0, WK model degenerates 
to ML model. Since fault-tolerant capacity is costly [oTfTUl 
[T7lfT8ll33] , a fundamental concern is how to efficiently allo- 
cate limited resources of capacity to makes network more 
robust. ML model raises a linear correlation between extra 
capacity and initial load, while WK model prefers to pro- 
tect the highest-load vertices. We agree with WK model 
that under the same cost, a certain more heterogeneous 
distribution of extra capacities can further enhance the 
robustness, however, the two-step function may be over- 
simplified. According to the previous works [10033] , the 
total cost e to protect a network can be denned as: 



1 N 



!)• 



(5) 



In ML model, the cost is 



N N 



eniL 



In the WK model, the cost is 



(G) 



i=0 



1 N 1 N T N" 



i=0 



i=0 



(7) 

where N" is the number of vertices with initial load larger 
than /?Lmax- 

In this paper, we propose a novel allocation mechanism 
of limited resources of capacity against cascading failures, 
as: 

Ci = (l + akl /{&))• I*, (8) 
where a > and 7 > are two free parameters, and 



1 N 



(9) 



The parameter a corresponds to the extra capacity, and 7 
is the parameter controlling the heterogeneity of resource 
allocation, the cost of this model is 



JV JV 



(10) 



Different from ML model, the fault-tolerant capacity 
and topological structure are interrelated in the current 
model. Since the components of a network may have dif- 
ferent fault-tolerant capacities, this work raises a more 
general allocation strategy for extra capacity. The assign- 
ing extra capacity depends not only on load Li of vertex i 
but also on the number of its links. When 7 > 0, vertices 
with larger degrees could be allocated more extra capac- 
ities. The case 7 = degenerates to ML model, and the 
case 7 — > 00 represents the extremely heterogeneous allo- 
cation where only the most connected vertex is protected. 
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Fig. 3. The robustness G vs. 7 under attacking of the highest- 
load vertex for UPTN on Beijing with size N = 4127 for dif- 
ferent a. 



Fig. 4. Go and Gmax vs. a. for UPTNs in Beijing (TV = 4127) 
and Shanghai (N=2035). 



3 Simulation 

Fig. 1 reports the robustness G as a function of 7 at sev- 
eral typical values of a for BA network [55114*0] of size 
N = 5000 and average degree (k) — 4. We find that 
the robustness has a maximum, G max , with 7 £ (0,1) 
for different a, indicating that the allocation of the lim- 
ited resource of fault-tolerant capability should be neither 
uniform nor extremely uneven, but a little bit more het- 
erogenous (for 7 > 0) than the case of ML model. When 
7 > 5, all G(j) at different a converge to a single curve, 
indicating that the extremely uneven distribution of re- 
source, with the largest-degree node possessing the ma- 
jority of extra capacity, has even worse performance then 
the uniform allocation (except for the cases with a close 
to zero). Define Go = G(j — 0), and Gmax = max 7 G(7) 
for given a. The former corresponds to the ML model, 
while the latter to the best robustness predicted by the 
current model. In Fig. 2, we plot Go and Gmax vs. a for 
the same data in Fig. 1. Clearly, in a broader region of a, 
Gmax is always greater than Go, demonstrating that we 
can get higher robustness of network than ML model with 
the same amount of resource. 

The serious global failure, as a catastrophe resulting 
from cascading propagation [13] . has often been observed 
in some kinds of traffic networks suffering overload failure, 
where the word "traffic" can stand for information pack- 
ets of the Internet [41]r42j , the electronic flow in the power 
grids [43], as well as passengers in public transportation 
systems [44,45] • Next, we compare the performance of the 
current model with ML model on some real networks. Es- 
pecially, the drawing near Olympic 2008 in Beijing and 
World Expo 2010 in Shanghai will bring extremely heavy 
traffic pressure on Beijing and Shanghai urban traffic sys- 
tems, thus an urgent problem is how to allocate the lim- 
ited additional traffic capacities to avoid possible cascad- 
ing traffic congestions. The topological characters of urban 
public traffic networks (UPTNs) are reported previously 



[46], and here we investigate the same data to see if our 
model can achieve better performance than ML model. 

In Fig. 3, we show the relation between the robustness 
G of Beijing's UPTN and the parameter 7 of allocation 
of capacity resources at various values of a. Qualitatively 
speaking, the pattern of G(7) is similar to that of BA 
network. Fig. 4 compared the performance of the current 
model and ML model for UPTNs in Beijing and Shanghai 
respectively. Clearly, the current model can always per- 
form better than ML model under the same amount of 
extra capacity. Moreover, large fluctuation is observed for 
ML model, thus one could not guarantee the enhancement 
of robustness giving more resource. In contrast, Gmax in- 
creases with a monotonously in the current model. Two 
G(a) curves in Beijing and Shanghai are remarkably dif- 
ferent for ML model, while almost have almost the same 
trend for the current model. Perhaps, this indicates that 
current model has some universality to improve robustness 
for different real-world networks, or maybe it is because 
the robustness of this strategy is near the optimal one. 



4 Summary 

In summary, the robustness of complex networks is ana- 
lyzed within the framework of allocating extra fault-tolerant 
capacity against cascading failure. We proposed a novel 
capacity model, which achieves higher robustness for the 
same amount of resource as ML model 10 . In previous 
works, the extra capacity of a vertex is linearly correlated 
with its load, or simply divided into two discrete values by 
a step function. In this paper, we argue that those prior 
strategies may be oversimplified, and a more complicated 
and more heterogenous (than ML model) allocation strat- 
egy could further enhance the robustness. The simulation 
results of B A networks and two typical real transportation 
networks strongly support the efficiency of our strategy. 
It is worthwhile to emphasize two features of the current 
model: Firstly, the maximum value of robustness Gmax 
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is a monotonously increasing function on parameter a, in 
contrast with ML model. Secondly, the same trend and 
pattern of G m ax(a) for different UPTNs indicates that 
this model may have some universality for robustness im- 
provement of different real-world networks (Maybe it is 
because the performance of this strategy is near the opti- 
mal one). However, this point is not clear thus far, which 
needs further investigation. Besides, although the previous 
works suggested using the mean relative extra capacity e 
to measure the cost, it is also, to some extent, reasonable 
to measure the cost by using absolute extra capacity. In 
that way, given a, the cost will be a monotonous function 
of 7, and our numerical results give raise to a practically 
significant conclusion: There exists a specific threshold of 
7 C corresponding to Gmax, below which one could improve 
the robustness by spending more, while above which the 
additional resource will, in contrast to what we expect, 
makes the network worse robust. Although at a first step 
far from the final goal, optimal allocation, we believe this 
model have its theoretical importance and potential appli- 
cation in designing infrastructure networks from the point 
of economic view. It can also provide guidance in designing 
more robust artificial networks. 
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